Density of states in d-wave superconductors disordered by extended impurities 
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The low-energy quasiparticle states of a disordered d-wave superconductor are investigated the- 
oretically. A class of such states, formed via tunneling between the Andreev bound states that 
are localized around extended impurities (and result from scattering between pair-potential lobes 
that differ in sign) is identified. Its (divergent) contribution to the total density of states is de- 
termined by taking advantage of connections with certain one-dimensional random tight-binding 
models. The states under discussion should be distinguished from those associated with nodes in 
the pair potential. 

PACS numbers: 74.25.-q, 74.72.-h, 74.25.Jb 



Introduction: In recent years, considerable attention has 
been focused on the low-energy electron-hole quasiparti- 
cle spectral properties of the cuprate superconductors in 
the presence of impurity scattering. Much of the impetus 
for this effort has its origin in the fact that many of the 
cuprate superconductors are randomly chemically doped 
insulators, and are therefore disordered. Moreover, as 
they are pair-breakers for them, the role of impurities 
is especially important for d-wave superconductors. Of 
particular interest is the behavior of the single-particle 
density of states (DOS) p(E) as the energy E tends to 
zero, i.e., its low-energy behavior. 

In recent work on the DOS of disordered d-wave super- 
conductors, Pepin and LeeEJ invoked a i-matrix approxi- 
mation to infer that p(E) ~ 1/Kllni? 2 | 2 at low energies. 
More recently, Yashenkin et al.El ancLAltlandB have ar- 
gued that the divergence found in Ref .EJ is present only for 
the case of a vanishing chemical potential (i.e. for a half- 
filled band), and thus does not apply to a doped cuprate. 
(To be precise, unitarity of the impurity .scattering is also 
required.) It is further argued in Refs£ra that, instead of 
diverging, p(E) should vanish at E = (unless certain 
very specific fine-tuning requirements are met). An im- 
portant feature shared by RefsBBB is the hypothesis that 
the disorder potential may be adequately modeled by a 
random collection of point-like scatterers. However, for 
a single impurity in a d-wave superconductor, the, low- 
energy DOS is qualitatively different for point-likco and 
extendedla (i.e. impurities of a size much larger than the 
Fermi wavelength) impurities: the states that reside at 
zero energy for extended impurities reside at nonzero en- 
ergies for point-like impurities. (The underlying reason 
for this difference is that for point-like impurities, the 
quasiparticle scattering is essentially diffractive, whereas 
for the extended impurities it is essentially semiclassi- 
cal.) This observation raises the possibility that such 
differences will continue to manifest themselves in the 
many-impurity setting. 

The purpose of the present Paper is to identify a mech- 



anism for producing low-energy quasiparticle states. This 
mechanism is based on impurity-scattering processes that 
connect states associated with differing signs of the 
wave pair potential^. In the case of a single impuritj 
this mechanism has already been shown to produce low- 
energy states that are localized near the impurity. These 
states can be associated with the classical trajectories 
scattering from the impurity, ancLhave been observed via 
scanning-tunneling spectroscopyuu. 

Here, we build upon this single-impurity physics to 
identify a singular (and potentially dominant) contribu- 
tion to the low-energy DOS. This contribution, which 
conventional techniques fail to capture, arises from tun- 
nelling along the classical trajectories that connect the 
individual impurities and, hence, connect the low-energy 
states localized near these impurities. The underlying 
physics was formulated some time ago in the context of 
tunnelling corrections to ground-state energips.in in mod- 
els of supersymmetric quantum mechanicsoiiljUj. 

The picture we have in mind(see Fig. Q) of the pro- 
cesses that lead to low-energy states involves classi- 
cal trajectories that each visit many extended impuri- 
ties. As a quasiparticle progresses along such a trajec- 
tory, its momentum is repeatedly altered via scatter- 
ing from the extended impurities, so that the effective 
pair-potential (Fig. |]b) undergoes sign changes. Local- 
ized near each such sign-change would be a zero-energy 
quasiparticle state; quasiparticle tunnelling (through the 
pair-potential) connects these states, raising their en- 
ergies from zero, and thus forming a_low-energy band 
that exhibits a Dyson-like singularity!! 2 ] at zero energy: 
p(E) ~ l/E\lnE 2 \ 3 . This picture loses its precision 
for sequences of impurities between which the pair po- 
tential is predominantly small (i.e. for nodal directions) 
and, consequently, the states are not well localized near 
the impurities. However, the contribution on which we 
are focusing (i.e. the non-nodal contribution) is expected 
to be substantial, and quite likely dominant, in the low- 
energy limit. 
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potential [i.e. k-p (= 27r/Ap) is the Fermi wave vector], 
V is the single-particle impurity potential, and we have 
adopted units in which h 2 /2m = 1, m being the common 
(effective) mass of the electrons and holes. The operator 
A is the pair-potential (integral) operator; how it acts 
is specified by the nonlocal pair-potential kernel A(r, r') 
via [A«„](r) = / dVA(r,r')u„(r'). 

To define this model fully, we need an assumption 
about the form of A(r, r'). It is convenient to exchange 
the coordinates r and r' for the relative and center-of- 
mass coordinates p and R: 

A(p,R) = A(r,r'), p = r - r', 2R = r + r' . (3) 



FIG. 1: (a) Sketch of a classical trajectory that encoun- 
ters several impurities and, hence, several sign-changes in A. 
The shaded regions on the trajectory denote the sections of 
the trajectory where approximate zero-energy states reside, 
(b) Momentum-space pair potential of a d-wave superconduc- 
tor, (c) Schematic depiction of A and A 2 along the trajectory. 

Based on previous work on low-energy quasiparticlc 
statcsoj, we identify three classes of processes that can 
effect the low-energy DOS of a disordered d-wave su- 
perconductor: (i) purely semiclassical scattering between 
states with differing signs of the pair-potential (i.e. scat- 
tering due to random extended scatterers in d-wave su- 
perconductors), which we shall focus on in the present 
Paper; (ii) purely diffractive scattering between states 
near the nodal directions of the pair-potential (i.e. the 
scattering of nodal quasiparticles by point-like_impu- 
rities), which was considered in, e.g., RefsBffiLj; and 
(iii) the mixing of the low-energy states that arise from 
processes (i) and (ii) . Throughout the rest of this Paper 
we shall ignore processes in class (iii). In consequence, 
the low-energy DOS of a disordered d-wave superconduc- 
tor can be expressed as a sum of contributions arising 
from class (i) processes (viz. p cxt ) and class (ii)r«rocesses 
(viz. pnodai), the latter, as discussed in Refs.Bu, being 
non-divergent and therefore subdominant. Thus, we shall 
focus on p cxt [and denote this by p{E)\. 
Quasiparticle motion in a d-wave superconductor with 
many extended impurities: Our focus will be on the DOS 
per unit area, i.e., 

where A is the area of the sample and the energy eigenval- 
ues {E m } follow from the Bogoliubov-de Gennes eigen- 
problem, viz., 

(ii) ft) =*»(:;)• < 2 > 

Here, h = — V 2 — k 2 , + ^(r), in which k 2 , is the chemical 



Then, via Fourier transformation with respect to p, viz., 



A(k, R) = J d 2 P e- lk "A(p,R), (4) 

we obtain the pair-potential A(k, R) at center-of-mass 
position R and relative momentum k. As our aim is to 
describe the cuprate superconductors, we take A(k, R) 
to have d-wave symmetry: A(k, R) oc (fc 2 — k 2 ,), where k x 
and k y denote the cartesian components of k. However, 
we allow for the possibility of position-dependent ampli- 
tude variations of the d-wave pair-potential due, say, to 
pair-breaking effects near the extended scatterers. 

Semiclassical approach to the Bogoliubov-de Gennes 
eigenproblem: We now invoke a semiclassical approxima- 
tion under which p(E) is expressed in terms of the solu- 
tion of a family of one-dimensional eigenproblems, each 
associated with a classical scattering trajectory in the 
presence of the single-particle impurity potential V(r). 
We restrict ourselves to a bejel, discussion of this ap- 
proach; for details see RefsMLj. The approximation 
amounts to our: (i) regarding the kinetic and potential 
energies as being comparable and being the largest en- 
ergies in the problem; (ii) turning off the pair potential; 

(iii) treating semiclassically the quasiparticle motion in 
the presence of the kinetic and potential energies; and 

(iv) reinstating the pair potential. Via this approach, we 
reduce the two-dimensional Bogoliubov-de Gennes eigen- 
problem to a family of one-dimensional Andreev eigen- 
problems residing on trajectories, each trajectory being a 
particular classical scattering trajectory in the presence 
of the many-impurity potential. This scheme applies un- 
der the following conditions: (i) the amplitude of A and 
V should vary slowly, relative to the Fermi wavelength 
Af; and (ii) the Fermi energy k\ should be large com- 
pared with the energy scale of interest, viz. E, as well as 
with the typical pair-potential scale. 

Let us now turn to the family of one-dimensional eigen- 
problems arising from this semiclassical scheme. Follow- 
ing Ref.B, these trajectory-dependent eigenproblems take 
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the form 



H 



E n 



_ f-2ik F 8 a A (ct) 
" ~ 1 Ao(a) 2ik F d a 



(5a) 

(5b) 
(5c) 



A (ct) ~ A(fc F <9 CT X c (cr),X c (cr)), 

i.e., an Andreev eigenproblemEl. Here, the parameter a 
measures the position along a particular classical trajec- 
tory x c (cr), the latter obeying Newton's equation in the 
many-impurity potential, viz., 



fcf ^xo(tr) = - W(x c (<r)). 



(6) 



The DOS is then obtained by assembling the eigenvalue 
spectra E n (n,b) of all the classical trajectories, the lat- 
ter being labelled in terms of an asymptotic momentum 
direction n and impact parameter b: 

P(E) * J dbp(n,b,E), (7a) 

p(n,b,E) = ^6(E-E m (n,b)). (7b) 

m 

Thus, in order to obtain p(E) one needs to find each 
classical trajectory, obtain the associated effective pair- 
potential [given by Eq. (Ml, solve the resulting one- 
dimensional eigenvalue equation and, finally, integrate 
over all the classical trajectories using Eq. ([7a]). We note 
that if one interprets the weight of a particular classical 
trajectory as the probability of finding a pair-potential 
configuration corresponding to that particular trajectory 
then we see that the calculation of p(E) amounts to com- 
puting the average density of states of a random pair- 
potential model. Models-Df this sort have been consid- 
ered, e.g., in RefsMBm. 

Eigenvalue problem for a single trajectory: We now ex- 
amine the contribution p(n, b, E) to the DOS for the case 
of a generic trajectory (n, b). For convenience, we intro- 
duce the rescaled trajectory parameter s = o/2k-p; the 
Hamiltonian then becomes 



H 



-id s A(s) 
A(s) id s 



(8) 



where A(s) = Ao(2fcps). 

Our method for calculating the spectrum of H in the 
many-impurity case is based on that for a single impu- 
ritycl. In the latter case, low-energy states arose from 
asymptotically sign-changing trajectories [i.e. those tra- 
jectories for which lims^ioo A(s) differ in sign]. Finding 
the spectrum amounted to identifying such sign-changing 
trajectories. What about the case of many extended im- 
purities? In this case, for a typical trajectory through the 
impurity potential A(s) undergoes repeated sign changes. 
On a particular trajectory let us label the the positions 
of these zeroes of A(s) by {s n }. Recall that we are con- 
cerned with the collection of impurity states that would 



lie at zero energy if the impurities were isolated. Ow- 
ing to tunnelling between them, these formerly degen- 
erate states yield a continuum of states, extending up- 
wards in energy from zero. Our task is to shed some 
light on this band-formation. We proceed to set up a 
tight-binding model along the trajectory, in which we re- 
tain only the zero-energy impurity states {|n)} (i.e. the 
local ground states at each of the {s„}) asd-allow only 
nearest-neighbor tunneling between themEJEil. To com- 
plete the model, we need the matrix elements of H con- 
necting these states, i.e. 



t n = (n\H\n + l) 



(9) 



Using thc_analytic expression for the zero energy wave- 
functionsEJ, (s\n) cx exp — ds'|A(s')|, a straightfor- 
ward calculation produces 

t n «-^|A'( s „)A'( S „ + i)| 1/4 exp{-^" +1 d s '|A( s ')|}, 

(10) 

where A'(s) = d s A(s). We are now in a position to write 
down a low-energy effective approximation to H, viz., 

£«^„(|n)(n + l| + |n + l)(n|), (11) 



i.e., for each classical trajectory, one arrives at a (topo- 
logically) one-dimensional hopping model that captures 
the physics of tunneling processes between the (formerly 
zero-energy) states localized near each zero of the pair- 
potential. 

Density of states: In order to obtain the low-energy DOS, 
we must obtain the DOS of the effective Hamiltonian jll| ) 
for each trajectory, and then collect them together. Wc 
assume that the collection of trajectories forms an ensem- 
ble that is characterized by the condition that momentum 
directions before and after a collision are uncorrelated. 
Then, summing over such an ensemble of trajectories is 
equivalent (up to a constant of proportionality) to aver- 
aging the DOS of the Hamiltonian ( pi] ) over uncorrelated 
values of t n . To obtain the low-energy DOS of this effec- 
tive modei-we appeal to results obtained by Eggarter and 
Riedingerca, who, building on the work of Dysont3 and 
Theodorou and CohenE3, studied random-hoppiag mod- 
els of precisely this form. Specifically, in RefB3 it was 
found that, under the condition that the {t n } are uncor- 
relatedcJ and identically distributed, the DOS as E — > 
is given by 



p(E)^N s Z 



2cj 2 



£|ln(£;/t") 2 | 3 ' 



(12) 



where N s denotes the average number of sites along the 
trajectory, Z is the constant of proportionality arising 
from the Jacobian of the transformation from summing 
over trajectories to averaging over {t n }, i is the scale 
characterizing {t n }, and the amplitude a 2 is given by 
the variance of the logarithm of t, i.e., 



<(ln* 2 /Ag) 2 ) 



(iniVA 2 ,) 2 . 



(13) 
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where (• ■ • ) denotes a disorder average. The scale t and 
a 2 could be estimated with the help of Eq. nG as 



t 2 cx k F A n l J 2 exp(-A /k F nl /2 ) 
a 2 cx A /k F nl /2 , 



(14a) 
(14b) 



where n c is the number of impurities per unit area. What 
remains is to determine the coefficient Z\ we now make 
an estimate of this quantity. 

Our estimate for Z follows from considering the inte- 
gration over b for a single impurity of size a. As the im- 
purity potential is expected to decay rapidly away from 
the impurity, the only trajectories that interact appre- 
ciably with it arc those that directly intersect it, i.e., 
k F J db — > k F a. Thus, in the absence of the tunneling 
corrections between the zero-energy states „p(E) is ap- 
proximately given by po(E) ~ k F an c S(E]B. Although 
the inclusion of tunneling corrections changes the energy 
dependence of p{E), here we assume that the energy- 
integrated density of states JJ p(E) dE is approximately 
conserved for some appropriately chosen cutoff e. These 
considerations lead to the following approximate form for 
p(E) valid for E < i: 



p(E)<x 



A an c 
E\\n(E/t) 2 \ 3 ' 



(15) 



We remark that the divergence of the low-energy DOS 
is ultimately cut-off due to physical processes not in- 
cluded in the present description. These include the de- 
phasing scale h/Tff, and the diffractive scattering scale 
S/ra.s. ; the cut-off will occur at the largest of these scales. 
Thus, it is possible that at extremely low energies p(E) 
eventually does vanish asymptotically, in agreement with 
the results of Refs.oEl. However, the present results would 
still apply at intermediate energies. Finally, we stress 
that the calculation presented here does not place any 
special emphasis on the nodes of the d-wave order pa- 
rameter; indeed, the singular contribution to the density 
of states of a d-wave superconductor reported here arises 
from states in generic (rather than nodal) regions on the 
Fermi surface. 
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